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abstract 

After a short review of the principal theories of spiral structure in galaxies, I describe two new devel- 
opments. First, it now seems clear that linear theory cannot yield a full description for the development 
of spiral patterns because V-body simulations suggest that non-linear effects manifest themselves at a 
relative overdensity of ^ 2%, which is well below the believed spiral amplitudes in galaxies. Second, 
I summarize the evidence that some stars in the solar neighborhood have been scattered at an inner 
Lindblad resonance. This evidence strongly supports a picture of spirals as a recurring cycle of transient 
instabilities, each caused by resonant scattering by a previous wave. 

Subject headings: galaxies: evolution - galaxies: haloes ~ galaxies: kinematics and dynamics - 
galaxies: spiral 



1. INTRODUCTION 

Despite many decades of effort, we do not yet have a 
widely accepted theory to account for the graceful spi- 
ral patterns in disk galaxies. Most workers in this field 
agree that spiral patterns are gravitationally driven varia- 
tions in surface density in the old stellar disk, as supported 
by photometric data (e.o. IZibetti et a/.ll2009l ) and stream- 
ing motions in hig h spatial resolution velocity maps {e.g. 
IShettv et al\\206lv . There is also general agreement that 
gas seems to be essential. 

There seems little doubt that some spiral patterns are 
tidally driven, while others could be the driven responses 
to bars. Although these two ideas may account for a large 
fraction of the cases, espec ially if orbiting dark matter 
clumps can excite patterns ()Dubinski et aL|[2008[ ). spirals 
can still develop in the absence of either trigger, as revealed 
in simulations. 

Here, I focus on the idea that spirals are self-excited 
oscillations of the stellar disk, which represents the great- 
est theoretical challenge. Two camps have long advance d 
quite distinct theories. One idea (e.g. iBertin fc Linlll996l) . 
is that spiral features are manifestations of quasi-steady 
global modes of the underlying disk. The other is that they 
are short-lived, recurrent transient pattern s that originat e 
either from shearing bits and pieces Je.g. 'To omrelll990l) . 
or something more subtle (e.g. ISellwo od 2000). 

A serious barrier to progress in this area has been the 
absence of observational discriminants that would favor 
one of these radically differing viewpoints over the other. 
The predictions for density variations or gas responses at a 
single instant are essentially independent of the generating 
mechanism and do not depend strongly on the lifetime of 
the pattern. 

In h is presentation, but not in the written version. [Toomrd 
()1981h likened theoretical work on spiral structure to the 
(apocryphal) blind men examining an elephant. Several 
important features of the beast have been described in pa- 
pers by many different authors over the past 50 years. Here 
I first review many of these components before attempting 
to weave them into a distinct, yet still incomplete, picture. 

The key idea described in this review resulted from a 
collaborative visit to Doug Lin in Santa Cruz in the sum- 
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mer of 1987. I developed it further for a few years, but 
worked on other problems whilst awaiting observational 
support for the somewhat complicated picture we pro- 
posed. The desired evidence appea red with the publication 
of the Geneva-Copenhagen Survey (jNordstrom et aZ.ll2004 
hereafter GCS), as I will explain below, and results from 
the Gaia satellite should improve matters still further. 

2. GLOBAL MODES OF ROTATIONALLY-SUPPORTED DISKS 

Simple mode ls of disk galaxies supp ort ma ny lin ear in- 
stabilities {e.g. iKorchagin et aLll2005t lJalalill200it) . The 

bar-forming mode is generally the fastest growing, but it 
has almost no spirality. These studies are therefore impor- 
tant to understand stability, but do not appear promising 
for spiral generation. 

The "den sity wave" theory fo r spiral modes, described 
in detail bv iBertin fc LinI (|1996D . invokes a more specific 
galaxy model with a cool outer disk and hot in ner disk. 
The local stability parameter, Q = (tr/ (iR^ciit (jToomrd 
11964 ). is postulated to be Q > 1 in the outer disk and to 
rise steeply to Q > 2 near the center. Under these spe- 
cific conditions, Bertin & Lin find slowly evolving spiral 
modes that grow by their WASER mechanism. They in- 
voke shocks in the gas to limit the amplitude of the slowly 
growi ng mode, leading t o a quasi-steady global spiral pat- 
tern. ILowe et al\ (I1994D present a model of this kind to 
account for the spiral structure of M81. The main objec- 
tion to their picture is that it is likely that such a lively 
outer disk will support other, more vigorous, collective re- 
sponses that will quickly alter the background state by 
heating the outer disk, as described below. 

3. RECURRENT TRANSIENTS 

V-body simulations of cool, shearing disks always ex- 
hibit recurrent transient spiral activity, and this basic re- 
sult has not change d for several d ecades as nume rical qual- 
ity has improved. ISellwood fc C arlberg ()1984D reported 
that patterns fade in simple simulations that do not in- 
clude the effects of gas dissipation; the reason is the disk 
becomes less responsive as random m otion rises due to par- 
ticle scattering by the spiral activity (jCarlberg &: SellwoodI 
[l985;.Binnev fc Laccv-1988.') . which is therefore self-limiting. 
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ISellwood fc CarlbereJ (|1984[ ) also show that mimicking the 
effects of dissipative infaU of gas, such as by adding fresh 
particles on circular orbits, allowed patterns t o recur "in- 
definitely." Later work (e.g. ^a rlber g fc FreedmanI Il985t 
lToomrdTl990l : IChakrabartil 120081) has shown that almost 
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any method of dynamical cooling can maintain spiral ac- 
tivit y, as also happens in galaxy formation simulations 
{e.g. iGovernato et aZ.ll2007[ ). 

4. INDIRECT EVIDENCE FOR TRANSIENT SPIRALS 

Thus the transient spiral picture offers a natural expla- 
nation for the absence of spiral patterns in SO disk galaxies 
that have little or no gas; maintenance of spiral activity 
requires a constant supply of new stars on near-circular 
orbits. There are several other indirect pieces of evidence 
that also favor the transient spiral picture. 

It has been clear for sometime that the velocity dis- 
persion of disk stars in the solar neighborhood rises with 
age (e .g. GCS) and also, for main sequence stars, with 
color (|Aumer fc BinnevI |2009[ ). which is a surrogate for 
mean age. The hi ghest veloc i ties cannot be produced 
by cloud scattering ( Lacevl[T99TI : iHanninen" fc Flvnnl[20Q2h 
and some other accelerati ng agent, such as trans ient spi- 
rals, seems to be required (|Binnev fc Lacevlll988[ ). 

Studies of the metallicit ies and ages of nearby stars 
(lEdvardsson et aLllT9 93: Nor dstrom eii a/.ll2004l : lReid et al\ 
l2007t iHolmberg et aLII2007r) find that older stars tend to 
have lower metallicities on average. As it is difficult to es- 
timate the ages of individual stars, the precise form of 
the relation is still disputed. However, there seems to 
be general agreement that there is a spread of metallic- 
ities at each age, which is also sup ported by other studies 
(jChen et al. 20031 ; iHavwoodl [20081 ). A metaUicity spread 
amongst coeval stars is inconsistent with a simple chem- 
ical evolution model in which the metallicity of the disk 
rises monotonically in each annular bin, without mixing 
in the radial direction (fSchonrich fc BinnevI l2009l) . Again 
ISellwood fc BinnevI (|2002t ) showed that the needed radial 
mixing arises naturally when the disk supports recurrent 
transient spirals. (See also Roskar et al. 2008ab). 

5. WAVE-PARTICLE SCATTERING 

Figure [I] shows, for an axisymmetric disk having a flat 
rotation curve, how the classical integrals of specific en- 
ergy, E, and angular momentum, L^, are changed for stars 
that are scattered by a steadily rotating mild potential per- 
turbation. The solid curve shows the locus of circular or- 
bits, which has slope flc, while stars with more energy for 
their angular momentum pursue non-circular orbits. The 
principal resonances for an m-fold rotationally symmetric 
disturbance that rotates at angular rate Vlp are marked. 
The familiar corotation and Lindblad resonances for near 
circular orbits can be generalized for non-circular orbits 
by using the angular frequencie s, fts fc Qr, of their dou- 
bly periodic motion defined in ([Binnev fc Tremainell2008l 
hereafter BT08, pp 146-147). For nearly circular orbits, 
flif, — >■ and flR — >■ K, which are the angular frequen- 
cies of circular and epicycle motion respectively. The dot- 
ted curves show the loci of the principal resonances where 
m(ri0 — Hp) — here Z = for corotation and I = ±1 
of the inner and outer Lindblad resonances respectively. 




Fig. 1. — The Lindblad diagram for a disk galaxy model. Circular 
orbits lie along the full-drawn curve and eccentric orbits fill the 
region above it. Angular momentum and energy exchanges between 
a steadily rotating disturbance and particles move them along lines 
of slope Qp as shown. The dotted and dashed lines show the loci of 
resonances for an m = 2 perturbation of arbitrary pattern speed. 



iLvnden-Bell fc KalnaisI (|l972l ) show that stars are scat- 
tered by spiral waves only at resonances. Since Jacobi's 
constant (BT08, eq. 3.112), is conserved in axes that rotate 
with the perturbation, the slope of all scattering vectors in 
this plot is AE/ALz — f^p. Since fic — Hp at corotation, 
scattering vectors are parallel to the circular orbit curve at 
this resonance, where angular momentum changes do not 
alter the energy of random motion. Outward transfer of 
angular momentum involving exchanges at the Lindblad 
resonances, on the other hand, does extract energy from 
the potential well that is converted to increased random 
motionj^ as is well known 



ISellwood fc BinnevI ()2002f ) show that scattering at coro- 
tation causes very effective mixing. In a few Gyr, multiple 
transient spirals cause stars to diffuse in radius over time. 
These changes at corotation of the spirals, which occur 
with no associated heating, are able to account for the 
apparent metallicity spread with age. 

Notice also from Fig. [1] that the direction of the scatter- 
ing vector closely follows the resonant locus (dotted curve) 
at the ILR only. Thus, when stars are scattered at this 
resonance, they stay on resonance as they gain random 
energy, allowing very strong scatterings to occur. The op- 
posite case arises at the outer Lindblad resonance (OLR), 
where the star quickly moves off resonance when it gains 
a small amount of angular momentum. 

6. SWING AMPLIFICATION 



The classic figure of dust to ashes presented bv iToonir^ 
(|l98lh and reproduced in BT08 (Fig. 6.19) shows the dra- 
matic transient trailing spiral that results from a small 
input leading disturbance. This linear calculation shows 
that the disturbance does not persist, but decays after its 
transient flourish. In the late st ages, the "wa ve action" 
propagates at the group velocity (|Toomrelll969l and BT08 
pp 499-500) away from corotation, where it is absorbed at 
t he Lindbl a d reso nances. 

iToomrel (|1990[ ) suggests that a large part of the spi- 
ral activity observed in disk galaxies is the collective re- 
sponse of the disk to shot noise in the density distribu- 
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Fig. 2. — The time evolution of the peak overdensity in a series of 
simulations of the half-mass Mestel disk with differ ent numbers of 
particles. The model is predicted by IToomrd l|1981f ) to be globally 
stable. The ordinate is the maximum values of 5 on grid rings. 
Linear theory predicts the amplitude should remain proportional to 
Ar-i/2. 

tion. iToomre fc KalnaisI ()1991[ ) show that ampHfied noise 
arising from the massive disk particles themselves can be 
understood in the shearing sheet, where the resulting spi- 
ral amplitudes are linearly proportional to the input level 
of shot noise. While this may be a mechanism for chaotic 
spirals in very gas rich disks, where the dumpiness of the 
gas distribution may make the seed noise amplitude un- 
usually high, it se ems likely that spiral amplitudes (e.g. 
IZibetti et a/.ll2009t ) are too large to be accounted for by 
this mechanism in most galaxies. Also, the spiral structure 
should be chaotic, with little in the way of clear symmetry 
expected. 

7. GROOVE MODES 

ISellwood fc KahnI ()1991| ) showed that disks are destabi- 
lized by a deficiency of stars over a narrow range of angular 
momentum, which creates a "groove" in a disk without 
random motion. The groove itself is unstable, and the 
instability becomes a global mode through the vigorous 
supporting response of the surrounding stellar disk. The 
resulting linear instability, which also develops in a disk 
with random motion, prod uces a large-scale spiral pattern. 
ISellwood &: BinnevI (2002) showed that the amplitude of 
the mode was limited by the onset of horseshoe orbits at 
c orotation. 

ISellwood k KahnI (|l99lD showed that almost any nar- 
row feature in the angular momentum density is desta- 
bilizing. Thus the common starting assumption of spiral 
structure studies, that the underlying disk is featureless 
and smooth, may throw the spiral baby out with the bath- 
water. 

8. A LINEARLY STABLE DISK 

The so-called "Mestel" disk has the scale-free surface 
density profile E — /{2nGR), with Vc being_the circu- 
lar orbital speed th at is independent of R. IZangj (jl976[ ) 
and iToomrd (|l98lh carried through a global, linear sta- 
bility analysis of this disk model with random motions 
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Fig. 3. — A power spectrum of m = 2 disturbances during the 
period of rapid rise in amplitude in the N = 5 X lO'^ experiment 
from Fig. [2] The solid curve shows 2Q{R) and the dashed curves 
2n{R)±K{R). 



described by a smooth distribution function. They in- 
troduced a central cutout and an outer taper in the ac- 
tive mass density, and replaced the removed mass by rigid 
components in order that the central attraction remained 
unchanged at all radii. The dominant linear instabilities 
they derived were confirmed in A'^-body simulations by 
ScUwood & Evans (20 0i|). 

iZanel (|1976D and lEvans fc ReadI (11998) showed that all 
full-mass disk models sufi^ered from global lop-sided insta- 
bilities no matter how large a degree of ra ndom mo t ion or 
however gentle the cutouts. Accordingly, IToomrd (|l98lt) 
halved the active disk mass and was able to prove that the 
resulting model was globally stable to all small-amplitude 
disturbances when Q = 1.5. This is the only known model 
for a non-uniformly rotating disk that is stable to small 
amplitude perturbations. 

Linear theory therefore predicts that A^-body simula- 
tions of the half-mass Mestel disk should exhibit no activ- 
ity in excess of the inevitable swing-amplified shot noise. 
Figure [5] shows that this is not the case. The ordinate is 
the largest value of (5 = AS/S from m = 2 disturbances 
in a sequence of simulations with increasing numbers of 
particles. The unit of time is Ri/Vc, where VcRi is the 
center of the inner angular momentum cut out. Thus the 
orbit period at this small radius is 27r. 

At t = 0, i5 oc N~^/^, as appropriate for shot noise, and 
swing-amplification causes an almost immediate jump by 
a factor < 10 for all N. When A^ = 5 OK, amplified noise 
causes ^ 20% overdensities, but the amplitude always rises 
to similar values in later evolution, when the inner disk has 
developed a pronounced bar, for all the particle numbers 
shown. A period of slow growth occurs in the simulations 
with the largest two values of A^, offering tentative support 
for the linear theory prediction of global stability, with 
the slow growth perhaps arising from the gradual growth 
of par ticle correlations as described bv lToomre fc KalnaisI 
(fT99l . Even in these cases, however, the amplitude rises 
more rapidly once i5 > 2%. 

Thus we always observe runaway growth of spiral ac- 
tivity - albeit more and more delayed as N is increased 
- behavior that is quite clearly not in accord with linear 
theory predictions. The rapid growth once 6 > 2% sug- 
gests that the behavior is already non-linear at this mod- 




Fig. 4. — The optimal coherent wave fitted to the m = 2 data 
from the simulation described in §1.5. The full-drawn circle marks 
the corotation resonance, the dotted circles the Lindblad resonances 
for this pattern. 



est amplitude. Note that the largest mimber of particles, 
N — 500M, is within a factor of 100 of the number of stars 
in a real galaxy disk, where in reality the mass distribution 
is far less smooth, owing to the existence of star clusters 
and giant gas clouds. 

9. MORE DETAILED ANALYSIS 

The runaway growth once S > 2% is of particular inter- 
est. A power spectrum analysis of non-axisymmetric den- 
sity variations during this period of rapid growth reveals 
multiple waves of many different frequencies, each extend- 
ing from the inner Lindblad resonance (ILR) to corotation, 
as shown in Fig. |31 We see that disturbances are present 
at several frequencies at any one time, and that each ap- 
pears and fades over a period of ^ 100 dynamical times, 
with those of lower frequency generally developing later. 

Figure 2] shows the transient wave that was prominent 
over a short time interval during which the wave grew and 
decayed while the pattern speed remained approximately 
constant. The measured fip — 0.182Vb/i?i, placing corota- 
tion and the Lindblad resonances for circular orbits at the 
radii shown by the circles in this Figure. The peak ampli- 
tude of the pattern is near corotation and the wave extends 
to the Lindblad resonances on either side; the outer half 
has lower amplitude because the wave is spread over a 
larger area. As the w ave decays, the corotation peak dis- 
perses fe ellwood fc B inncv 2002) and the "wave action" is 
carried radially at the group velocity where it is absorbed 
at the Lindblad resonances. 

It is convenient to separate the specific energy of a star 
in an axisymmetric potential, z), into three distinct 
terms E = Ec\rc{L^) + E'ran + E^, where 

Scire = $(i?c,0)+L2/(2i?2) 

i?ra„ = $(i?,0)-$(i?„0) + i vl + Ll(^^^ 
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Fig. 5. — The distribution of particles in the space of angular 
momentum and energy of non-circular motion at t = 600 in the 
simulation mentioned in the text. The resonant loci (solid) and 
scattering trajectories (dashed) are also shown. 



Here Rc is the radius of a circular orbit having the same 
Lz, and vj^ is the radial velocity. The approximate form 
for Ez (BT08, eq. 3.86) holds for stars that do not climb 
to great heights above the mid-plane, with v being the 
vertical oscillation frequency of such stars. 

The change in random energy caused by an angular mo- 
mentum change AL^, is AE'ran = ^^pAL^ — EdidLz + 
ALz) + Ecirc{Lz). It is easy to show that scattering from 
a circular orbit according to this formula takes place along 
a line in (L^, £'ian)-space that is has a steep negative slope 
for the ILR and positive slope for the OLR. 

Figure [S] shows the distribution of particles in the space 
of angular momentum (here labeled J) and energy of ran- 
dom motion (here labeled E — Ec) in this simulation at 
t — 600. It differs from the initial distribution by the 
pronounced tongues of particles at low angular momen- 
tum that reach up to high random energies on an upward 
curving, negatively sloped trajectory that is almost com- 
pletely obscured by the loci of the Lindblad resonances 
(solid lines), and the scattering trajectories (dashed lines). 
These lines are for the measured pattern speed of the wave 
shown in Fig.|4]and therefore have no free parameters. The 
fact that the dashed line at the ILR overlays the tongue 
of particles so perfectly that they are barely visible indi- 
cates that these particles were scattered to more eccentric 
orbits by the ILR of the identified pattern. Recall from 
Fig. [T] that particles stay on resonance as they gain angu- 
lar momentum, which is the reason some reach very high 
energies. 

10. RECURRENT CYCLE? 

Each coherent wave leaves behind an altered distribu- 
tion function and apparently thereby creates the condi- 
tions fo r a new instability. De spite having shown this al- 
ready in lSellwood fc LinI (|1989[ ). I am still unclear how this 
happens in detail. I have not put much effort into pursuing 
this question since then, because I was concerned that the 
behavior I was studying might turn out to be a horrible 
artifact of the simulations and thought it prudent to await 
evidence that something similar occurs in nature. I hoped 
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Fig. 6. — The distribution of GCS stars in the space of the two 
integrals -Eran and L^. Aside from the general decrease in density 
towards large -Eran and a skew to lower Lz caused in part by the 
asymmetric drift, the DF in this projection also shows significant 
substructure. 



()Sellwo od 1994") to see evidence of resonance scattering in 
the HIPPARCOS measurements of the local stellar kine- 
matics. 

11. GENEVA-COPENHAGEN SURVEY 



iNordstrom et al\ (j2004l ) followed up the accurate dis- 
tances and proper motions from HIPPARCOS to obtain the 
missing radial velocities from ground-based spectra. Their 
survey has yielded known distances and full space motions 
of over 14000 local F & G dwarfs - a sample that is be- 
lieved to be almost complete within 40 pc and largely free 
from selectio n biases to 100 pc. I u se the latest revision of 
their survey (jHolmberg et a/.ll2009l) . discarding only stars 
having distances > 200 pc and all known Hyades cluster 
stars, which leaves a sample of 13 045 stars. 

The addition of radial velocities confirmed the extraor- 
dinary substruc t ure in the U — V plane already discov- 
ered bv iDehnenI ()1998[ ) in the absence of radial velocities. 
The local distribution function is far from the smooth dou- 
ble Gaussian postulated by Schwarzschild (BT08, p. 323); 
there is little evidence for an underlying smooth compo- 
nent, rather the distr ibution has the a ppearance of sev- 
eral distinct streams ("Bow et al. 2009), which can not be 
dissolved clust ers (Bensby et al. 2007; F amaev et al. 2007: 
iBovv fc Hog'3[20090 . Various groups have attempted to 
account for this substructure by resonant s cattering by the 
Milky Way bar (Dchncn 2000), spir al arms (iDe Simone et al\ 
[200l . or both ( Quillen fc Minchevl[2005H Ghakraba rtvl2007D . 
which may indeed account for parts of the observed struc- 
ture. But none of these studies was self-consistent, and all 
attempted to explain the structure in velocity space. 

It is more instructive to project the data into integral 
space, as shown in Figure IHl The distribution in the up- 
per panel assumes a locally flat rotation curve for the 
Milky Way, i?o = 8 kpc and Vq = 220 km/s, but the 
features in the plot are insensitive to these assumptions. 
The almost parabolic lower boundary is a selection effect, 
because stars whose angular momenta differ increasingly 
from that of the LSR do not visit the Solar neighborhood 



unless they have greater energies of random motionQ The 
left-right bias results from the fact that the density of stars 
in the Galaxy declines outwards, and therefore larger num- 
bers visit the solar neighborhood from the innner Galaxy 
than from larger radii {i.e. the asymmetric drift). 

The distribution in Figure [5] manifests one strong tongue 
with a steep negative slope rising from an abscissa Lz — 
0.96i?oVb. This c ould be the scat tering feature at an ILR 
that I suggested (jSellwoodl |1994[ ) might be present, but 
in fact trapping at any of the principal resonances would 
cause an excess of stars along a line of similar slope. Nev- 
ertheless the feature is a highly significant density ex- 
cess, with ^ 9 9% confidence from a bootstrap analysis 
(|Sellwoodll20Tol) . While there are hints of other features, 
none but the obvious one is statistically significant. 

In order to determine which resonance could have caused 
this feature, we must examine the phases of the stars, 
which should exhibit different correlations depending upon 
the resonance responsible. I t herefore compu ted action- 
angle coordinates for each star (|Sellwoocill2010D , and show 
the results in Figure [71 The radial coordinate in the upper 
panel shows the approximate epicycle size a = \/2 J/j/k 
while the azimuth is 2w^ — wr. The concentration of stars 
at one phase, marked by the arrow is highly significant. 
Furthermore, the stars colored red in the skinny triangle 
in the upper panel are exactly the stars, also colored red in 
the lower panel, that defined the scattering tongue. This 
Figure therefore confirms that the tongue of stars in inte- 
gral space has been caused by scattering at an ILR. 

Unfortunately, it is not possible to tell from this analy- 
sis the angular periodicity of the disturbance that created 
this feature; it could be an m = 2 spiral, but higher an- 
gular periodicities would give rise to very similar features 
in these plots. However, an m = 2 wave seems the most 
likely, since the alignment of the scattering trajectory and 
resonance locus in Fig. [T] is closest for this angular peri- 
odicity, making it easier for strong scattering to produce 
a clear feature. 

The identification of an ILR provides strong support for 
the picture of spiral generation that I have been developing 
from the simulations. It suggests that spirals in the Milky 
Way are transient, with the decay of one pattern seed- 
ing a new instability tha t produces another s piral pattern. 
An ILR also excludes the lBertin fc LinI ()l996l ) quasi-steady 
wave hypothesis as the cause of the pattern that gave rise 
to this particular feature, although this one case clearly 
does not preclude the possibility that it may still apply 
elsewhere. 

12. CONCLUSIONS 

Globally stable A^-body simulations have manifested re- 
current transient spiral patterns for many years, and the 
phenomenon has not changed as numerical quality has 
risen. The transient nature of spiral patterns receives indi- 
rect support both from the importance of gas and secular 
heating of disks. The more recently discovered further 
conse quence of radial mixing in disks ()Sellwood fc BinnevI 
I2OO2I) is important for structural evolution, metallicity gra- 

^ A similar excluded region is visible in the quantities plotted by 
Arifyantom & Fuchs ( 2006); the omission of the small correction for 
the different Galacto-centric distances of the guiding centers in their 
energy-like term is unimportant. 
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Fig. 7. — Top panel shows the distribution of GCS stars in action- 
angle space. The radial coordinate is the square root of the radial 
action, while the azimuthal coordinate is 2w^ —wn. Note the excess 
of stars with phase ~ 53° (marked by the arrow). The lower panel 
shows the distribution in (X/z,£'ran) space. The red points are the 
same stars in each plot. 



dients, dynamo theory, etc. 

Simulations of the only known stable disk with non- 
uniform rotation have revealed that linear perturbation 
theory breaks down at relative overdensities of just ^ 2%. 
Real galaxy disks, which contain both star clusters and 
molecular clouds, are not that smooth and it would seem 
therefore that their behavior cannot be fully described by 
linear theory. I show how non-linear effects could lead to 
a recurrent cycle of spirals, although the exact mechanism 
remains obscure. 

The Geneva-Copenhagen survey has full phase space co- 
ordinates of an unbiased sample of local F & G dwarf stars. 
Examining the distribution as a function of integrals has 
yielded clear evidence for the existence of an inner Lind- 
blad resonance in the solar neighborhood. This discovery 
provides strong support for the mechanism for recurrent 
spiral activity that I have been developing. 
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